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We calculate the contribution of deconfined spinons to inelastic neutron scattering (INS) in the 
fractionalized antiferromagnet (AF*), introduced elsewhere. We find that the presence of free spin- 
1/2 charge-less excitations leads to a continuum INS signal above the Neel gap. This signal is found 
above and in addition to the usual spin-1 magnon signal, which to lowest order is the same as in the 
more conventional confined antiferromagnet. We calculate the relative weights of these two signals 
and find that the spinons contribute to the longitudinal response, where the magnon signal is absent 
to lowest order. Possible higher-order effects of interactions between magnons and spinons in the 
AF* phase are also discussed. 



I 

O 

o 



> 

00 

m 
o 

O 



o 
o 



X 



I. INTRODUCTION 

Theories of spin-charge separation in the high-T^ 
cuprates have been hotly debated almost since the orig- 
inal discovery of these materials [Q. Finding a the- 
ory of electrons in more than one spatial dimension 
which exhibits zero-temperature spin-charge separation 
has proved to be as theoretically challenging as it is phe- 
nomenologically appealing 1^. It can be argued that all 
such theories will admit, in the low-energy limit, a formu- 
lation in terms of a Z2 gauge theory ||3|. Recent papers 
have addressed the problem of finding microscopic mod- 
els of electrons which become fractionalized in some range 
of their parameters Q. It remains an important task to 
enumerate concrete, experimentally- measurable conse- 
quences of these exciting theoretical ideas. Previously, 
we have explored the consequences of two-dimensional 
fractionalization on the spectral function, as probed by 
angle-resolved photo-emission spectroscopy In this 
paper, we calculate the inelastic neutron scattering signal 
from spinons in a fractionalized antiferromagnet (AF*). 
We find that these spin-1/2, charge-less excitations lead 
to a continuum of excitations above a gap. Because we 
are interested in the parent insulators of cuprate super- 
conductors, we have taken a phenomenological model 
for the spinons which gives them both a Neel gap aris- 
ing from antiferromagnetic ordering and a rf-wave pair- 
ing gap which becomes the pseudogap at moderate dop- 
ing and the superconducting gap in the superconduct- 
ing phase. We contrast this signal with the signal from 
excitations in a conventional antiferromagnet and calcu- 
late the strength of the spinon signal compared to the 
magnon signal (which is also present). This comparison 
estimates the feasibility of measuring this anomalous sig- 
nal in the parent insulators. We also discuss higher-order 
effects stemming from interactions between spinons and 
magnons. 



II. THE MODEL 

The AF* phase has been discussed elsewhere |^-||] 
and here we use the same phenomenological model in- 
troduced and justified |^,^ previously. We assume that 
the steps of: (1) deriving a lattice Hamiltonian contain- 
ing all important effective interactions between electrons 
and (2) splitting the electron into chargon and spinon 
fields (cia = bifia) and deriving the appropriate Z2 gauge 
theory have been performed and we have arrived at the 
following effective low-energy Hamiltonian for the system 
in 2-dimensional fractionalized phases: 
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where the spinon pairing Ay 

+A along X 
—A along y. 
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is taken to be d-wave: 



Ay = 



(1) 
(2) 

(3) 



and the spin operator is S,; ~ ^fjcrfi- Here, <i,j> are 
nearest neighbors on a 2d square lattice. The U term is 
a Hubbard- like interaction for (I ~ x) chargons per unit 
cell. 

We now briefly justify this model for the underdoped 
cuprate materials on phenomenological grounds. For suf- 
ficiently small doping and low temperatures such that the 
Z2 theory exhibits fractionalization, the Hamiltonian is 
as written in Eq.(0). At temperatures below the energy 
scale A, the spinons arc effectively paired into rf-wave 
singlets and there is a d-wave gap to spin-1/2 excita- 
tions. For large enough g (and an additional minuscule 
3d spin coupling) the system develops long-range anti- 
ferromagnetic order. At half-filling, the chargons are 
gapped into an insulating phase and we obtain a frac- 
tionalized insulator with long-range Neel order and an 
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additional d-wave gap to spin-1/2 excitations, previously 
dubbed AF* Moving away from half-filling, the an- 
tiferromagnetic order will be quickly suppressed, while 
for tc U and with an additional long-range Coulomb 
interaction, one still expects the chargons to be insulat- 
ing. We then have a fractionalized insulating phase with 
a d-wave gap to spin-1/2 excitations. Within a spin- 
charge separation scenario, this phase is identified with 
the pseudogap regime in the cuprates. For chargon hop- 
ping, tc, sufficiently large, the chargons Bose condense, 
giving a d-wave superconductor. At large dopings, we 
expect the system to recover Fermi liquid properties, as 
occurs when the vortex excitations (visons) in the Ising 
gauge field condense thereby confining the spinons and 
chargons to form the electron. A schematic phase dia- 
gram is shown in Fig. ^ In this paper, we elucidate 
further some of the properties of the AF* phase, found 
at half-filling. 

Recent experiments by Bonn, Moler,ef al put limits 
on the likelihood of this sort of spin-charge separation 
in Y Ba2Cu'iOQ+x, although the experiments have only 
been performed on one sample so far The question 
of spin response in an antiferromagnet which is fraction- 
alized is nevertheless well-posed and could be relevant to 
other materials. Also, it is quite possible that some other 
sort of exotic order lurks in the cuprates; this work would 
then serve as an illustrative calculation. 
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FIG. 1. Schematic phase diagram for the high Tc cuprates 
within a spin-charge separation scenario. 



III. EFFECTIVE HAMILTONIAN FOR THE SPIN 
SECTOR 

In this paper we work at half-filling, where the charge 
degrees of freedom will be gapped into a Mott insulating 
phase, and calculate the spin response of the system, ap- 
propriate for magnetic probes such as neutron scattering. 
Hence, from here on, we assume that the relevant piece of 
the Hamiltonian in Eq. is that containing the spin de- 
grees of freedom and we ignore the charge degrees of free- 
dom. At temperatures much less than the vison energy. 



the chargons and spinons are essentially non-interacting, 
so this is reasonable in a fractionalized phase. 

Hg (Eq.(^)) may be decoupled in a path integral, us- 
ing a Hubbard-Stratonovich transformation. This gives 
us the following low-energy theory for the spin sector: 

Hsp.n = i-t-flfl + Au/a/^i + H.C.] 

<i,j> 

-g J2 N,:, p • (S, - S,;+^) + I ^ (N,,^)^ (4) 

where N is a 3-component vector of classical fields liv- 
ing on the nearest-neighbor links of the lattice, which we 
have broken into its two square sublattices, labeled A and 
B and shown in Fig. ^. G {±x, ±y}. 
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FIG. 2. The 2d square lattice with sublattices A and B 
marked. 

We begin by analyzing Eq.(^ at the mean field level. 
We find antiferromagnetic Neel ordering and a quadratic 
theory for the spinons which can be solved exactly. We 
derive a self-consistency equation for the mean field 
approximation which connects the magnitude of the 
Neel order to properties of the spinons. Fluctuations 
about this mean field solution lead to spin-waves (as in 
the conventional antiferromagnet). Thus, we find that 
the fractionalized antiferromagnetic phase has two spin- 
carrying excitations: the spin-1 magnons and the spin- 
1/2 spinons, which interact with each other. To deal 
with these interactions, we treat the spinons as a per- 
turbation: setting ts — A = 0, one recovers the con- 
ventional antiferromagnet and its spin- wave excitations. 
Integrating out the spinons to each order in tg/g and A/g 
would give modifications of the spin-wave theory due to 
the spinons. If one wishes to find higher-order proper- 
ties of the spinons, one may start with the mean field 
result and then integrate out the magnons, generating 
interactions between the spinons. In the limit is, A <C g. 
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it is clear that the theory in Eq.(Q) can be solved in a 
controlled fashion. 



A. Mean Field Theory 

First, we ignore the spinons, setting ^ A — 
and reducing Si to the usual spin-1/2 quantum opera- 
tor. Eq.(0) becomes: 



'9 



(5) 



Choosing z as the spin quantization axis, this is mini- 
mized classically by: 



Si I ^)ieA - 2I ^)^eA, 



(6) 
(7) 



leading to a Hamiltonian for the field N: 

^f-fEE(N.^.)'-5EEN.-.- (8) 



i^A fi 



i£A M 



This is minimized by: (N|^) = Nq = 1 yi,fi. 

Plugging this mean field solution for N back into 
Eq. gives a Hamiltonian for the spin sector at the mean 
field level which is quadratic in the spinons: 



H 



MFT 



<i,j> 



-4g^iVoi-Sr + | E ^0, (9) 

i i£A,fj, 



where we have retained z as the spin quantization 
axis and have written = X^i ^7 = 

5 X^k fw+Tv^f^- lattice spacing has been set to unity. 

The full solution to this quadratic spinon Hamiltonian 
has been given elsewhere [ppO|] and here we reproduce 
only the dispersion: 



with 



Ng = 2gNo, 

ek — — 2ts (cos kx + cos ky ) , 
Ak = — A(cos kx — cos ky). 



(10) 



(11) 
(12) 
(13) 



B. Self-Consistency of the Mean Field Solution 

The mean field solution with A^o = 1 is found for- 
mally in the limit tg — A — Q, and one expects the 
spinons to reduce the Neel order from this maximum 
value. We therefore take a mean field solution of the 
form (N| ) — Nq and demand that it be a saddle point 



of the full theory, Eq. 



Z[No] 
Z[No] 



'[No] 



-0H'^^'^[NoJ,r] 



EE e"'^ 



(14) 
(15) 



Tracing out the quadratic spinons gives an expression for 
5"^[iVo] and the saddle point condition, SS'^'^/SNq = 0, 
gives a self-consistent equation for A^o- At zero tempera- 
ture, this is: 



1 = 2.9 



d^k 



1 



(16) 



(27r)2 £k[NoV 

with Sk[No] given in Eq.(^. 

It is easy to check that for tg = A = Q, Eq.(p^ gives 
A^'o = 1. For ts, A ^ g, one finds the perturbative result: 



Nn 



, {2tsr+A^ 



(17) 



At large values of tg/g and A/g, one expects the spinons 
to drive the Neel order to zero, even at zero temperature. 
In order to keep our calculations controlled, we work in 
the limit ts,A <C g and treat the spinons as a pertur- 
bation. The reasonableness of this limit for the physical 
systems in question will be discussed later. 



C. Fluctuations About the Mean Field 

We see from the self-consistent mean field calculation 
of the previous section that the Neel order is reduced at 
zero temperature by the spinons. Even in the absence of 
the spinons (i.e. in the pure-spin model with tg = A = 0) 
we know that fluctuations in the order parameter are im- 
portant and lead to magnons. This suggests the follow- 
ing program for calculating the spin excitations in AF* 
beyond the mean field level. First, set = A = and 
work with the spin-1/2 quantum operator, S^. "Integrat- 
ing out" these operators on each site leads to an effective 
theory of fluctuations in the field N and gives the spin- 
wave dispersion. Then, one may integrate out the spinons 
perturbatively in tg/g and A/g. This will lead to inter- 
actions between the magnons and give corrections to the 
magnon dispersion. 

With tg and A set to zero, the effective spin Hamilto- 
nian, Eq.(^), reduces to Hf[N, S], Eq.(^). In this section 
we look at fluctuations of the field N around the mean 
field solution. We therefore set N^ = 1 on all links and 
look at the fluctuations, N — zN^ ~ N^. Hg'^ can then 
be written in the following form: 
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ieA M 
i£A jeB 



(18) 

(19) 
(20) 

(21) 



Integrating out the operators S,; on each site amounts to 
performing perturbation theory in iJ'"'. To second order, 
the resulting effective action for N-'- is: 



(22) 



i£A,^i 



i ieA 



(e EAN.!.P + E EANi--.,.p) 



(23) 



where r is imaginary time and we have set the tempera- 
ture to zero. 

Since we are interested in obtaining a long-wavelength 
theory for the spin-waves, we make the coarse-graining: 



N^, i £ A. This amounts to working in the 



basis of the Goldstone modes of the theory. With this 
approximation, we arrive at the effective action: 



+ E 



E 



dT 



16<? 



, (24) 



where all sites are on the A sublattice, {■■■) refers to 
nearest-neighbor pairs, and ((■■■)) refers to next-nearest- 
neighbor pairs. Fourier transforming to k and (imagi- 
nary) Lo gives: 



S'S= / 5|N^(k,c.)| 



^1 + 7^) 



with 



7k = -j^icosk^ + cos ky). 



(25) 



(26) 



This immediately gives the lowest-order magnon dis- 
persion, LUk = ^gy/ (1 — 7fe)/(l + 7fe)- This dispersion 



is similar to the usual one found using, for instance, 
Holstein-Primakov bosons (see, for instance, Ref. p|), 
but the ratio of the spin-wave velocity to the maximum 
of the dispersion is different. Even neglecting spinons en- 
tirely, this is only a lowest-order result for fluctuations of 
the field N. Working to higher orders in the perturbation 
theory of Eqns.(p^pl|) will generate a more realistic spin- 
wave theory of the conventional antiferromagnet. For the 
purposes of calculating the INS response, we satisfy our- 
selves with the lowest-order result for the magnons. 
The lowest-order effect of the spinons, as we have seen 



in Section IllB, will be to reduce (N^) from unity. At 



higher order, we expect interactions with the spinons to 
further affect the magnon dispersion. However, it would 
be wise for the purposes of comparing the magnon and 
spinon INS responses to take into account this over-all 
reduction of the Neel order by the spinons. Therefore, 
we make the substitution: g gN^, where A^o is set 
by the self-consistency condition, Eq.([l^), giving a "self- 
consistent" magnon dispersion: 




(27) 



One important consequence of this calculation is 
that it allows one to set the parameter g in terms of 
experimentally-measured properties of the magnon spec- 
trum. Because it is probably less sensitive to details of 
magnon interactions, we choose to set g by the maximum 
of the magnon dispersion rather than by the spin-wave 
velocity near k = (0,0). In the Heisenberg model, the 
maximum of the magnon dispersion is 2 J; we therefore 
make contact with this model by identifying J — gNo. 
Many experimental probes of the undoped cuprates find 
J to be around 150 meV 



IV. INELASTIC NEUTRON SCATTERING 
RESPONSE 

What is the difference in spin response between the 
model given in Eq.(^) and a spin-charge confined an- 
tiferromagnet? The lowest energy spin excitations in 
both systems are the spin-1 magnons dictated by Gold- 
stone's theorem in this symmetry-broken phase. For 
the confined insulator, the lowest-energy spin-1/2 exci- 
tations would presumably be something like single elec- 
trons, which have a huge Mott gap (on the order of a 
few eV in the cuprates). In contrast, we see immedi- 
ately that provided one is in the regime where Eq.(0) 
holds (i.e. at temperatures and energies small compared 
to the vison gap), the lowest energy spin-1/2 excitations 
are spinons which propagate as independent excitations 
above the Neel (and rf-wave) gap which is on the order of 
J ~ 0.1 eV. In the previous sections, we have laid out a 
theory of the spin degrees of freedom in a fractionalized 
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insulator with long-range Neel order. Now, we calculate 
the INS signal in AF* using the lowest-order theories of 
the magnons, Eq.pS), and spinous, Eq.(ra). 



A. Magnetic Neutron Scattering Cross-Section 

The differential cross-section for neutrons scattering by 
wave- vector q = k/ — k.^ and energy oj off electronic spins 
is Hi 



a, 13 



(a,/3 — X, y, z), where i^^(q) is a form factor and the spin 
structure factor is: 

s'-'i^M = ^ ,_^i/...T i"^x"^(q--)- (29) 



universal features. By universal properties we mean that 
the limits of both the response function and the spin- wave 
dispersion as q ^ (0, 0) and q (tt, tt) are the same for 
all calculational methods because they are dictated by 
symmetries. Including intra-spin-wave interactions mod- 
ifies the non- universal aspects of the dispersion {e.g. the 
spin-wave velocity near q = 0) and yields direct multi- 
magnon contributions to the spin structure factor. These 
direct multi-magnon processes are of a much lower weight 
than the single magnon processes and so we ignore them. 

To first order then, the magnon response to neutron 
scattering in AF* has the same universal features as in 
a conventional antiferromagnct. This calculation allows 
us to fix the parameters in our theory: we have seen in 
Section tllC that the energy scale gN^ — J to this or- 
der. Additionally, the magnitude of the direct spinon 
signal can be compared with the magnitude of this well- 
established magnon signal. 



where 

X"'^(q,zc.„)= rrfre'""^(T,S^(r)S^q(0)) (30) 
Jo 

is the imaginary-time spin-spin correlation function. At 
temperatures such that fc^T -C we can take the zero- 
temperature limit: 



1 



I _ ^-Lj/ktT 



ei-Lo) 



(31) 



with 8(x) the Heavyside-step function. The spin-1/2 
operators are given by the usual expression with elec- 
tron operators replaced by spinon operators: S(q, r) = 

The model outlined in the sections above provides a 
theory of the spin response in AF*. In the next sec- 
tion, we use the lowest-order results to calculate the 
magnon and spinon signals (respectively) in inelastic neu- 
tron scattering. We include a brief discussion of higher- 
order effects. 



C. Spinon Response 



As detailed in Section III A, at the mean field level 



the spinon part of the Hamiltonian in Eq. (|1|) is quadratic 
and has been solved elsewhere The spin-flip and 

longitudinal structure factors for the spinous are as fol- 
lows: 



1 - 



a 



1 



fq-k^^k 

: 57+(q,w), 



x,5(c^-£q_k-fk), (33) 



1 - 



Eq-k 



(34) 

Aq_kAk 
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^q— k / \ <^k / ^q— k^k 



-l-elastic (Bragg) response. 



(35) 



B. Magnon Response 

Starting from the ts — A ^ spin Hamiltonian, 
Eq.(||), it is straightforward to calculate the spin-spin 
response function by including a source term in the ef- 
fective action of the form: ' '^i- The result, to 
lowest order in fluctuations of N is: 



^magnons 



(l-7q)^ 



(1 + 7^)^^134 



S{LO-io^), (32) 



with 7q and uiq given in Eqns.(|26|-p7D. While the exact 
form of this response function differs from that found in, 
say, the Holstein-Primakov formalism, it has the same 



where Ng, efc, A^, and £k are defined in Eqns. ( [10[]13| ). 

These formulas have a few salient features. Spin-fiip 
neutron scattering leaves a spin-1 excitation in the sam- 
ple with momentum q and energy uj. The expression 
for iSy " simply sums the ways of destroying a spin-down 
(-up) spinon at momentum — k and creating a spin-up 
(-down) spinon at momentum q — k, with the constraint 
that the energy cost for this process must be lo. The rest 
of the expression is the zero-temperature probability that 
the state at — k is occupied and the one at q— k is unoccu- 
pied, appropriate for fermions. If one takes Ng = (cor- 
responding to no Neel order), the full spin-rotation in- 
variance of the spinon system is restored and 5+" = S^" , 
as expected. 
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For INS using unpolarized neutrons, the differential 
cross section is a combination of tfiese two signals, but 
it is possible using polarized neutrons to obtain signals 
from these two channels separately (albeit at a large 
cost to the intensity). It is worth noting that, pro- 
vided one aligns the polarized neutron spins along the 
direction of the staggered magnetization, the signal due 
to single magnons is only in the "spin-flip" sector, i.e. 
'^f magnon " 0. The first contributiou to S^^ from spin- 
waves alone occurs in the two-magnon channel which 
is of a greatly reduced weight compared to the single 
magnon signal. The spinous then constitute the only 
magnetic signal of reasonable weight in this channel at 
energy scales of order J. 




N_0 



FIG. 3. Results of numerically integrating the RHS of 
Eq.(|l|), called I{No), for t^/g = A/g = 0.7. The 
self-consistent value of No is the one for which I — 1. 



The above expression for 5y ~ can be numerically in- 
tegrated to obtain the spinon response. This has been 
performed, approximating the (5-function in energy with 
a Lorentzian of width e = O.OIJ. None of the salient fea- 
tures of the results were influenced by the speciflc values 
of parameters (such as f^. A, e). 



A discussion of parameters is wise at this point. Work- 
ing in units of gNo — J, we will use tg = A = J to cal- 
culate S"^ for the spinons. This gives numerical values 
of these constants which are reasonable for the undoped 
cuprates. However, one may wonder whether they violate 
our above working assumptions ts/g,A/g ^ 1. Indeed, 
it is not immediately obvious how to set the parameter 
g. Here, we take the following tack: the quantity gN^ 
is set by comparison between the magnon dispersion in 
Eq. ( P7|) and the maximum of the magnon dispersion 
in INS experiments. This gives gNo — J ~ I'aQmeV . 
The parameters ts and A may similarly be set by exper- 
iments and we use the values given above. With these 
two constraints, the self-consistency equation, Eq. (16), 



becomes an equation for the parameter g. It is clear that 
if the values of ts and A are too large, the equation for g 
will have no solutions. For the parameters above, using 
the expansion (valid for small ts/g and A/g) in Eq.(|l^) 
as a first iteration yields g ~ 1.4 J. Plugging the resul- 
tant values of ts/g and A/g into the full self-consistent 
equation for N^^ the right hand side of Eq.(|l6|) can be 
numerically integrated for various values of Nq. The re- 
sult is shown in Fig. |^. It is clear from the graph that 
the value of A'o which satisfies the self-consistent equa- 
tion is A^o — 0.7. This gives gN^ = .98J ~ J, and we 
have our self-consistent parameters. We note that the 
value A'o — 0.7 is still rather close to the "no spinon" 
mean field value of unity, so that treating the spinons as 
a perturbation is somewhat justified WA- 



D. INS Signal 

In Figs. P and ^ we present the lowest-order results for 
— S^^ + '^liagnona J corrcspondlug to spin-fiip neutron 
scattering. All energies are in units of J = gN^. Note 
that we have taken the staggered magnetization to point 
in the z-direction, but since our theory does not contain 
terms which couple the spin and spatial variables (as, 
say, a spin-orbit coupling would), the spin axes can be 
rotated independently of the spatial axes. 



ji/2 




FIG. 4. Contours in the q^, qy plane. 
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As 5"^(q, is a function of three variables (for the 
effectively two-dimensional cuprates of interest), we show 
here a few views of this function. Fig.^, shows contour 
plots of the intensity as a function of distance along cuts 
in the {qx,qy) plane (shown in Fig.^ and energy. For 
the magnon signal, the 5- function in Eq.(p2[) has been 
replaced with a "box" function: 



Six) 



1/e for -e/2 <x < e/2, 
else. 



(36) 



with e = 0.06J ~ lOmeV. In these plots, we see the 
magnon dispersion with zeros at q = (0, 0) and (tt, tt) 
and a vanishing weight near (0,0). Above this, we see 
the spinon continuum turn on with a lower-bound which 
is modulated with twice the period of the magnon dis- 
persion. We see in this plot that the spinon signal is 
appreciable. 



(a) 







FIG. 5. 5+~(q,t^) along cuts (a),(b),and (c) m Fig. |. En- 
ergies are in units of J. White regions are zero intensity and 
black regions are high intensity. 



(5E.0) 



In Fig.|^, we show contour plots of the spinon intensity 
as a function of and qy at constant values of the energy. 
We see that the spinon intensity at turn-on [lo = 4J) is 
largest near the corners at (0, 0), (tt, tt), etc. For compar- 
ison, Fig. ^(d) shows plots of the single magnon weight 
(Eq.(^2[)) as a function of q^ and qy. 




(a) 
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FIG. 6. Spinon S^^{q,ui) at energies oj/J = (a) 4.5 (b)5.0 
(c)6.0 . The single magnon weight, Eq.(|32|), is shown in (d). 
White regions are low intensity and black regions are high 
intensity. 



In Fig. 1^ we present contour plots of S'|^(q, to) sequen- 
tially along the cuts in Fig. ^. We note again that, to 
first order, there is no magnon contribution to this corre- 
lation function and so the two-spinon continuum visible 
in these plots should in principle be the primary magnetic 
source of INS in this channel. 



S-- {ci^ a>) 




2 
1 

[0,0) tlC.Oj (H^lCj [3t/2,lLaJ (7E.0) 

FIG. 7. 5"(q,ij) along a cut {a.)+{h) + {c) in Fig. §. En- 
ergies are in units of J and the color scheme is the same as in 
Figs. I and ^ 



E. Higher Order Effects 

At higher orders in the perturbation theory of Eqns. 
( |l9| , we would expect spin wave interactions (present 
and important for detailed features even when tg and A 
are zero) to modify the magnon dispersion and lead to 
multi-magnon signals in both the 5+" and S^^ channels. 
This calculation would be the same as for a conventional 
antiferromagnet (since is = A = 0) and would presum- 
ably lead to the actual spin-wave response in a conven- 
tional antiferromagnetic system. At higher orders in tg/g 
and A/g, we would obtain magnon interactions mediated 
by spinous, which are not present in conventional anti- 
ferromagnets. Being gapped excitations, the spinous are 
inherently a "high energy" phenomena as far as the spin- 
waves are concerned. We might therefore expect them 
to influence most heavily the high-energy portions of the 
magnon dispersion. One can see from the graphs given in 



Section IV D that the lower edge of the spinon continuum 
has a minimum near (tt, 0), where the magnon dispersion 
has one of its maxima. One would expect interactions be- 
tween these two excitations to effect the dispersions most 
heavily near these points of closest approach. While the 



effect of unconventional excitations such as spinons can 
in principle be seen in the spin-wave response, we have 
chosen in this paper to focus on the direct spinon contri- 
bution to inelastic neutron scattering, shown in the pre- 
vious section. The spin-waves will also affect the spinon 
signal at higher orders, however, probably not in any 
dramatic fashion. 



V. CONCLUSIONS 

In this paper we have presented a theory of the spin 
excitations of the fractionalized antiferromagnet, AF* , in 
the limit where the visons can be ignored. This theory 
is Eq.(^). This phase has both long-ranged antiferro- 
magnetic order and the topological order associated with 
fractionalization [ p^ . It contains two spin-carrying exci- 
tations: the spin-1 magnons and the spin-1/2 spinons, 
which interact with each other. The theory is well- 
defined and can be solved in a controlled manner in the 
limit ts/g,A/g <C 1. In this paper, we have found the 
lowest-order theories of these two excitations and have 
calculated the dynamic spin-spin response functions of 
each, appropriate for inelastic neutron scattering experi- 
ments. 

To lowest order, the magnon signal is the same as in 
a conventional antiferromagnet, but higher-order effects 
of the spinons should lead to modifications of the dis- 
persion, particularly near (7r,0). The main anomalous 
feature of the INS signal from AF* is the presence of a 
spinon continuum at energies ~ 4J, which exists even 
in the longitudinal response, where the magnons are ex- 
pected to be absent at lowest-order. 
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